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We adapt to neutrino physics a general formulation for particle propagation in fluctuating 
media, initially developed for applications to electromagnetism and neutron optics. In 
leading approximation this formalism leads to the usual MSW effective hamiltonian gov- 
erning neutrino propagation through a medium. Next-to-leading contributions describe 
deviations from this description, which arise due to neutrino interactions with fluctuations 
in the medium. We compute these corrections for two types of fluctuations: (i) microscopic 
thermal fluctuations, and (ii) macroscopic fluctuations in the medium's density. While the 
first of these reproduces standard estimates, which are negligible for applications to solar 
neutrinos, we find the second can be quite large, since it grows in size with the correla- 
tion length of the fluctuation. We consider two models in some detail. For fluctuations 
whose correlations are extend only over a local region in space of length I, appreciable 
effects for MSW oscillations arise if (8n/n) 2 £ > 100 m or so. Alternatively, a crude model 
of helioseismic p- waves gives appreciable effects only when (Sn/n) > 1%. In general the 
dominant effect is to diminish the quality of the resonance, making the suppression of the 
7 Be neutrinos a good experimental probe of fluctuations deep within the sun. Fluctua- 
tions can also provide a new mechanism for reducing the solar neutrino flux, giving an 
energy-independent suppression factor of \, away from the resonant region, even for small 
vacuum mixing angles. 



1 



1. Introduction and Summary 



We do not understand our sun as well as we should. Although the sun shines brightly 
in neutrinos, it is not so bright as it ought to be according to our present understanding 
of its workings and of neutrino properties. 

In recent years experimental and theoretical lines of research have converged to bring 
this Solar Neutrino Problem (SNP) to a head. On the one hand, confidence in the 
experimentally-measured neutrino fluxes has grown with the observation of the solar- 
neutrino shortfall in four independent experiments [1], including those which are capable 
of detecting neutrinos from the principal p — p cycle of nuclear reactions. Moreover, it 
has recently become possible to calibrate these detectors by exposing them to very intense 
radioactive sources here on earth. On the theoretical side, confidence in the neutrino-flux 
predictions of solar models has also improved, for two reasons [2]. First, the redundancy 
of the experiments permits the discrepancy to be mainly based on the predictions for the 
p — p neutrinos. Since the p — p reactions are largely responsible for generating the sun's 
energy, the theoretical uncertainty in their reaction rates is minimal. Second, the rise of 
the field of helioseismology has made available an abundance of experimental data about 
the solar interior, thereby significantly improving the constraints on the assumptions which 
must be made in constructing solar models. 

If the problem is not the sun, then the measured neutrino shortfall must arise while 
the neutrinos are en route to the earth. Besides gaining support from the improvements 
in understanding of solar models, the credibility of such a neutrino solution to the SNP 
is also boosted by the existence of a very plausible and elegant mechanism for depleting 
the observed solar neutrino flux. The mechanism consists of resonant (MSW) oscillations 
of the neutrinos as they pass through the sun [3]. In this picture the small influence of 
the solar medium on neutrino propagation plays an important role by resonating with 
the equally small vacuum oscillations which generically arise once neutrinos are endowed 
with masses. Considerable effort has been invested in understanding the nature of these 
material-dependent oscillations. 

A common feature of the majority of these studies has been the approximation in 
which the influence of the solar medium is described in terms of an effective hamiltonian, 
depending on the mean values of the quantities to which the neutrinos couple. Less has 
been done to study scattering from the deviations away from this mean. Some researchers 
have investigated the effects of neutrino scattering from position-dependent densities [4], 
although usually ignoring the potentially decohering effects [5] - more about which later 
- of such scattering. Incoherent scattering due to interactions with the particles which 
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make up the medium has also been studied within the context of supernovae, and the 
early universe, [6] [7] [8] , in which case neutrinos themselves can be part of the underlying 
medium. This type of scattering is entirely negligible within the sun [9]. Until quite 
recently, [10] [11] [12], less attention has been devoted to the effects for solar neutrinos of 
more macroscopic fluctuations. 

The purpose of the present work is to develop a framework for describing the influence 
of all such fluctuations on neutrino propagation, with the goal of identifying when each 
can be important. To this end we adapt to neutrino physics a formalism which has been 
extensively used to describe the interaction of electromagnetic waves and neutrons with 
fluctuations in matter [13]. As we describe in detail herein, our results agree with earlier 
approaches when applied to the fluctuations they consider. 

In order of magnitude, fluctuation effects contribute to neutrino evolution with strength 
G\ {8nSn) Here the average is over different quantities in different situations, and 
5n = n — (n) denotes the deviation of the particle density from its mean. £\\ is the cor- 
relation length along the direction of neutrino motion. The relative size of this term as 
compared to the usual MSW evolution term, G F (n), is therefore of order G F {n) £\\e 2 , 
where (5n5n) = e 2 (n) 2 . Using the central solar density, (n) ~ 10 26 /cm 3 , we see that 
sizable effects can be expected only for large-scale fluctuations: £\\e 2 > (G F (n)) -1 ~ 100 
km. Our more detailed analysis shows that for resonant oscillations this estimate is too 
large, and sizable effects can arise starting from e 2 £\\ > 100 m. These scales are of potential 
interest for solar neutrinos, since they are typical of scales which can arise from physics 
within the sun. 

Some work on this kind of macroscopic-fluctuation-driven effects for neutrino propa- 
gation has appeared recently in the literature, starting with the pioneering work of ref. [5] . 
This reference considered a fixed density profile which varied in space, and computed the 
time evolution of the reduced density matrix which governs the flavour degrees of freedom. 
Decohering fluctuation effects were found when the neutrino momenta were integrated out 
to obtain the flavour evolution. We argue in Section 5 that this type of fluctuation is 
unlikely to be important for solar-neutrino physics, although it could well play a role in 
other applications. We reach this conclusion because we find that the only fluctuations 
which can decohere neutrinos as they pass through, are those whose size, £±, transverse to 
the direction of neutrino propagation, is smaller than the transverse size of the detector. 
But since fluctuations in the sun can in any case only affect neutrino evolution for £\\ > 100 
m, correlation lengths as small as typical neutrino detector sizes can play no significant 
role, so long as £± ~ £\\. The same need not be true for other applications, such as to 
supernovae, however. 
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More recently, the main approach to fluctuations that has been pursued to date, e.g. 
by refs. [10] and [12], is to model the density as a Gaussian random variable subject 
to the assumption that all correlations exist only over distances that are negligibly small 
compared to the neutrino oscillation lengths which are of interest. Our results extend these 
analyses in several ways. First, since we work from first principles, we can give explicit 
expressions which are applicable to any kind of density ensemble. In particular, we do 
not assume the correlation length of the fluctuations to be small compared to neutrino 
oscillation lengths, and so can apply our results to density profiles which vary on scales 
that are comparable to the size of the sun. We can also incorporate arbitrary variation in 
space and time of the fluctuation's mean and variance. This permits us to consider such 
real density variations as helioseismic waves, which are known to exist in the sun. Our 
equations reduce to those of refs. [10] and [12] in the limit that our assumptions overlap, 
but some of our most interesting applications are to situations for which previous analyses 
do not apply. 

Our results are presented in the following way: 

• 1: We present a general formalism in Section 2 for describing the interactions of any 
particle with arbitrary matter fluctuations. One of the main features of such fluctuations 
is that they generically destroy the coherence of neutrino propagation by evolving pure 
states into mixed states. As a result it is typically impossible to describe them in terms 
of a matter-dependent effective hamiltonian, since any such hamiltonian would necessarily 
take pure states to pure states. This section culminates in a master formula for the rate of 
change of the density matrix describing particle propagation in an arbitrary medium, which 
naturally divides into a term which defines a mean effective matter-dependent hamiltonian, 
plus a fluctuation-dependent term. A nice feature of the formalism is its recursive nature, 
which permits fluctuations to be successively integrated out on larger and larger distance 
scales. 

• 2: Section 3 applies the general results of Section 2 to neutrinos moving in the presence 
of microscopic fluctuations, which are those for which the correlation length is negligible 
compared to the interesting scales for neutrino propagation. We consider in some detail 
the special case of thermal fluctuations, and rederive the usual result that these are small 
for neutrinos in the sun. 

• 3: Section 4 then considers matter fluctuations on larger scales. These fluctuations arise 
because, although any one neutrino sees a fixed density profile, successive neutrinos see 
different ones. We describe this by considering neutrinos to pass through an ensemble of 
density profiles. The ensemble is characterized by expanding the density in terms of a 
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complete set of modes whose amplitudes are taken to be uncorrelated random variables. 
The nature of the underlying physics governs the basis of modes which are uncorrelated 
in any given application. We consider two types of bases for illustrative purposes: (i) 
fluctuations which are localized in position, within cells of slowly-varying length, £; and 
(ii) fluctuations in the amplitude of the normal modes which describe acoustic density 
waves. In the limit of small, constant, I the first of these bases reduces to the case studied 
in refs. [10] and [12]. The second is new, and is meant as a crude model of a helioseismic 
p-wave of fluctuating amplitude. We evaluate, for both examples, the mean hamiltonian 
and the contributions of fluctuations to neutrino evolution. 

• 4 : Section 5 takes the previous results and integrates out the neutrino momentum 
degrees of freedom, to obtain the reduced evolution equation which governs the reduced 
density matrix describing neutrino flavour and spin. We argue that previous derivations 
[5] , which trace over neutrino momenta without taking into account that neutrino positions 
are ultimately measured in real experiments, give mistakenly large estimates of the size of 
the decoherence which this trace introduces for solar neutrinos. A general expression is 
found for the fluctuation contributions to the neutrino-flavour evolution, which are found 
to be characterized by a single parameter, A a b- This parameter is evaluated for the two 
density ensembles introduced in Section 4. 

• 5: Section 6 specializes the general results to the two flavour case, and integrates the 
time evolution to obtain the electron-neutrino survival probability. An approximate an- 
alytic form for this integration is obtained, which is a generalization of Parke's formula 
for standard MSW mixing. The decoherence due to fluctuations appears in the evolution 
as a damping term, similar (but not identical) to what would happen if the neutrinos 
were decaying. Numerical integration is also performed, and found to agree well with the 
analytical results. 

We perform the MSW analysis for both models of density fluctuations that are de- 
scribed in Section 4. It is found that appreciable changes to the usual MSW scenario 
arise for surprisingly small amplitude fluctuations. In terms of e 2 = (5n5n) / (n) 2 , and 
the correlation length £y, we find deviations for e 2 £y > 10 m. Startlingly large changes 
from MSW behavior arise for e 2 £y > 1 km. A generic new feature of fluctuations is the 
introduction of a universal energy-independent reduction of the survival probability to | 
for small E/5m 2 . 

The model of a fluctuating helioseismic p-wave gives discernible effects which are 
comparatively small. For a wave with a 30 minute period, discernible effects require 
e > 1%. We understand the size of this effect to be due to the small wavelength the 
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wave typically has near the resonance region, due to the increase of the speed of sound 
with depth in the medium. More realistic simulations are currently under way to see if the 
same is true for neutrinos propagating through both solar p and g waves. 

• 6: Finally, in Section 7 the general formalism is applied to derive the effective evolution 
equation governing the reduced neutrino flavour /spin density matrix in the presence of a 
magnetic moment interaction. Our conclusions are briefly summarized in Section 8. 

2. The General Formalism 

In order to keep all approximations explicit it is instructive to first formulate our 
problem within its most general context. Suppose, therefore, that our system consists of 
two sectors, A and B, of which we wish to follow the evolution of degrees of freedom in 
sector A while ignoring (or partially ignoring - see below) those in sector B. For example, 
when examining the influence of matter fluctuations on neutrino propagation we will take 
A to describe the neutrino states of the system while B consists of the states which are 
available to the electrons and/or nucleons which make up the medium through which the 
neutrinos move. 

It is sometimes necessary to consider the more general case where a partial measure- 
ment is made on sector B, in addition to the measurements which are performed in sector 
A. We do this in order to set up the treatment of resonant oscillations, for which A con- 
sists only of the flavour (and spin) sectors of the single-particle neutrino sector, while B 
contains both the neutrino position/momentum information as well as all medium-related 
effects. (We argue in Section 5 that an improper treatment of this case has in the past led 
to a mistaken estimate of the size of incoherent effects purely due to this removal of mo- 
mentum degrees of freedom.) The slightly more general formulation is required to analyze 
this situation since neutrino position information is in practice never completely ignored 
(i.e. neutrinos are all detected on Earth). 

At an initial time, t', we suppose these two sectors to be completely uncorrelated. 
That is, suppose the initial density matrix for the entire system factorizes: 

p(t') ~~ g a ® Q B - (1) 

We imagine here that g A acts only in the A sector of the Hilbert space and g B acts only in 
the B sector, and we take g A and g B to be separately normalized within their own sectors: 
Tr A g A = Tr s g B = 1. Here Tr A (or Tr B ) denotes a trace taken only over the A (or B) 
sector of the Hilbert space. 
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Next suppose the hamiltonian for the system takes the form: 



H = H + V, 



(2) 



in which H = H A + H B describes the separate evolution of the A and B sectors, while 
V is the interaction which couples these sectors together. Using this hamiltonian we may 
evolve p to later times, t > t'. We assume for this purpose that g A and g B respectively 
commute with H A and H B . Within the interaction representation the time evolution of 
p(t) is then described by: 

dp 



dt 



— i 



(3) 



where V(t) = e iH()t V e~ iHot . Alternatively: 



with 



p(t) = U(t,t')p(t')U*(t,t') 
dU(t,t') 



dt 



= -iV(t)U(t,t'). 



(4) 



In general this time evolution will introduce correlations between sectors A and B and so 
won't preserve the factorized form of eq. (1). The remainder of this section is devoted to 
explicitly displaying these, and other, effects as sector A evolves in the presence of sector 
B. 

2.1) Coherent and Diffuse Scattering 

Suppose, now, that only observables associated with sector A are to be measured at 
some time t > t' . The probability of the results of any such measurement are completely 
described by the reduced density matrix, p A (t), defined by tracing the full density matrix 
over only the B sector of states: 



P ^)eeTt p(t) 



(5) 



Notice that eq. (1) implies p A (t) satisfies the initial condition p A (t = t') = q a . 



7 



We now wish to split the time evolution for p A (t) into a piece which describes the 
mean features of sector B — 'coherent' scattering 1 - - plus a piece which describes the 
fluctuations about this mean — 'diffuse' scattering. Our guiding principle in so doing is 
to ensure that final time-evolved probabilities may be written as the non-interfering sum 
of a coherent part plus a diffuse part. 

Define, then, the mean (or coherent) evolution operator, U(t,t'), as the average of 
U(t,t') over the B sector, as follows: 



U(t,t') = (U(t,t')} 



B ' 



(6) 



where the S-average of any quantity is defined by: (• • ■) B = Tr B g B (- ■ •) 
between U(t, t') and U(t, t') we denote: 



The difference 



AU(t,t') = U(t,t') - U(t,t') 



(7) 



and so satisfies the defining identity (AU) B = 0. 

With these definitions all probabilities calculated at times t > t' are the sum of a 
coherent piece and a diffuse piece. That is, for any hermitian observable acting only in the 
A sector, A , eqs. (1), (4), (6) and (7) imply: 



(O a ) (t) = Tr p(t)0 A 
= Tr 



U(t,t')p{t')U*(t,t')0 A 

U{t,t')Q A U*{t,t')0 A ] +Tr 
= (0 A ) c (t) + (0 A ) d (t). 



= Tr 

A 



AU(t,t')p{t')AU*{t,t')0, 



(8) 



The cross terms involving both U(t, t') and AU(t, t') vanish by virtue of the identity 
(AU) B = 0. This last equality defines the mean (or coherent) and fluctuation (or diffuse) 
parts of (O a ) (t). 



1 We borrow the descriptions 'coherent' and 'diffuse' from the analogous applications of this formalism 
to the propagation of X-rays and neutrons through matter. This split, as made precise in eq. (8), is our 
definition of coherent scattering for the present purposes. Notice that the 'coherent' part, as defined 
here, is coherent only in a weaker - though more useful, for present purposes — sense than is sometimes 
used in electromagnetic applications. For instance, coherence in the present context need not imply phase 
coherence between incident and scattered waves. 
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The distinction between diffuse and coherent evolution can also be made directly for 
the reduced density matrix itself. That is, p A (t) = p c A {t) + p d A {t), where 



p A (t) = U(t,t')g A U*(t,t'), 



P d At) 



Tr 

B 



AU{t,t')p(t')AU*{t,t') 



(9) 



For many applications — including the description of neutrino oscillations — it is 
preferable to formulate the diffuse-coherent split for dp/dt rather than for the integrated 
evolution operator, U(t, t'). This is because it is often possible to use perturbation theory 
for dp/dt but not for the long-time evolution of p(t). This leads us to formulate the main 
result of this section. The differential evolution equation for p c A {t) may be written: 



9Pa 
dt 



V(t)p A (t) - p A (t)V (t) 



(10) 



where V(t) is the effective interaction hamiltonian which is defined in such a way as to 
ensure that V is related to U(t, t') in the same way that V is related to U. That is: 



in) 



which need not be hermitian (since U need not be unitary). 
Similarly, the differential evolution equation for p A {t) is 



dpi 

dt 



d m 
dt b 



AU(t, t')p(t')AU* (t, t')J 



(12) 



2.2) Incorporating a Partial Measurement in Sector B 

For some applications it is true that measurements do not completely ignore what is 
going on in sector B. For instance, in applications to solar-neutrino oscillations we will 
follow common practice and take A to describe only the neutrino flavour and spin degrees 
of freedom. This involves banishing all neutrino position and momentum information into 
sector £?, even though any realistic measurements do include some information concerning 
neutrino position, such as that they are detected on earth. This section describes the slight 
generalization of the formalism which is required to handle such cases. 
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We therefore relax the assumption that all of the observables of interest, O, need act 
only in sector A. Instead we assume them to involve a specific observation in sector £?, 
which is uncorrelated with all measurements in sector A. That is, consider the class of 
observables having the form: 



o = o A <g> o. 



(13) 



with O a (or O b ) acting only in sector A (or B). In this case expressions similar to those 
found above may be derived, in which all averages over sector B are weighted by the 
observable O b . 

Specifically, define once more the evolution operators U(t,t') and AU(t,t') as in 
eqs. (6) and (7), but with the S-average now defined by: 2 



(-) B = 





'{■■■)q b O b 


Tr B 


q b O b 





(14) 



This choice preserves the property that (O) (t) may be written as the non-interfering sum 
of a diffuse and coherent contribution, although eq. (8) is slightly modified to become: 



(O) (t) = Tr p(t)0 



Tr 



U(t,t')g A U\t,t')0 A 



Tr 

B 



QbO, 



+ Tr 



AU(t,t')p(t')AU*(t,t')0 



(15) 



(0) c (t) + (0) d (t). 



As before, the time evolution of any such observable may be completely described in 
terms of a reduced density matrix, p A (t), for which the definition, eq. (5), is now replaced 
by: 



p A (t) = Tr p(t)0 



(16) 



satisfying the initial condition: p A (t = t') = q a Tr b [q b O b ]. Notice that p A (t) defined this 
way is not normalized. The differential time evolution of its coherent part is now given by 
the analog of eq. (10): 



dt 



V(t)p A (t)-p c A (t)V*(t) 



Tr 

B 



QbO l 



(17) 



Beware: the operator ordering in this definition has the counterintuitive implication that (X)* B need 
not equal (X*) B . 
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Notice that this equation lacks a term proportional to ^§f-, even though O s (t) generally 
depends on time within the interaction picture with which we are working. Its omission 
from eq. (17) is justified since there O b should be evaluated at the time, t m , when the 
measurement is performed, rather than at the time, t, of the evolution. V(t) is once again 
defined by eq. (11). 

The diffuse evolution now becomes: 

AU{t,t')p{t')AU*{t,t')0 B ^. (18) 



op 



at 



u m 
T^Tr 

8t b 



2.3) Perturbative Expressions 

It is instructive to evaluate eqs. (11) and (18) perturbatively in the interaction V. To 
this end we use the familiar series solution to eq. (4): 



* ft fT n -l 

U(t, t') = VH) n / dn • • • / dr n V(n) ■ ■ ■ V(r n ). 
n=o J * Jt ' 



(19) 



Using this expression in the previous results leads to the following formula for effective 
hamiltonian: 



V(t) = (V(t)) B -ifdr (SV(t) SV(t)) b + O (V 3 ) , 
Jt' 



(20) 



with SV(t) = V(t) - (V(t)) B . 

Notice that the antihermitian part of V first arises at second order in V. For instance, 
if [ e „ O,] = 0: 



1(V + V')=(V) 

W- r ) = -\[ dT {{mt),m.T)}) B + o(v>) 



-t 



5V(t),5V(r) ) +0{V 3 ), 



(21) 



Similarly, the rate of change of p d A is: 



dp d 



-q± = J'dr ^[{5V{t)p{t')5V(r) + 5V{r)p{t')5V{t))0 B ] + O (V 3 ) . (22) 



Eqs. (20) and (22) are our starting point for applications of this formalism to neutrino 
propagation through matter. 
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2.4) Long-Time Evolution and Master Equations 



Regardless of how small the interaction hamiltonian should be, perturbation theory 
eventually fails if one follows the system's evolution for sufficiently long times. Worse, it is 
often precisely the long-time behaviour which is of interest in particular applications. In 
this section we outline how to use the above perturbative expressions for time scales which 
are sufficiently large compared to the correlation times which govern the fluctuations in 
the medium. 

Suppose, then, that the correlation, {5V(t)5V(t')) : is negligible for \t — t'\ greater than 
some correlation time, r. Suppose also that the system's time evolution is required over 
timescales, T, for which T»r. In this case perturbative expressions for ^f- may be useful 
provided that r is small enough to justify perturbation theory, even if the same would not 
be true for timescales as large as T. 

In this limit a coarse-grained time derivative of the density matrix may be defined for 
times which are large compared to r [13]. It is given by neglecting the difference between 
p A and p c A in the expressions for -j^, as well as neglecting the difference between p(t) and 

f) rfi 

p{t') in Since these differences are higher order in the perturbation, V, this neglect is 
justified over time scales which are short enough to lay within the domain of perturbation 
theory. With these approximations, the perturbative expression for the sum of eqs. (10) 
and (12) becomes: 



dp A = dp c A dp d A 
dt dt 



dt 

v 2 (t)p A (t)-p A (t)v;(t) 



(23) 



+ 



dr Tr 

B 



6V(t) p(t) 6V(t) + 8V{t) p(t) 6V(t))o B + O (V 3 ) , 



where V 2 (t) denotes the second-order expression, eq. (20). 

If the correlation scale, r, is now assumed to be small compared to the time scale over 
which the coarse graining is taken, then we may neglect correlations on the right-hand-side 
of eq. (23), by writing p(t) « p A (t) <g> p B . With this choice eq. (23) describes a Markov-like 
process, for which ^f- depends only on p A (t), and not on the behaviour of p A for times 
previous to t. This represents a great simplification once eq. (23) is integrated to obtain 
the evolution of p A for very long times. It is this form of the time-evolution equations 
which is used in Section 6 for describing neutrino evolution within the sun. 
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2.5) Unitarity and Decoherence 



There are two general features of the above expressions which bear special emphasis. 
Notice first that the condition Trp = 1 implies the same is true for the reduced density 
matrix (when O b = I): Tr A p A = 1. This is easily seen to follow from eqs. (10), (20) and 
(22) by virtue of the following identity, which expresses the optical theorem in the present 
example: 



^p A = -t^{p A (t) [V(t) - V*(tj\ } + Tr 



dpi 
dt 



= 0. 



(24) 



Second, dpi/dt and V — V both cause a loss of coherence within sector A. That is, if 
the system is initially prepared in a pure state, for which p A = p A , then it need not remain 
so after interacting with sector B. The endpoint of evolution is therefore generally a mixed 
state. Quantitatively, starting from an initially pure state eqs. (10) and (12) imply the 
following rate of coherence loss: 



d ( 2 



P 2 ,=PA ^ J 



P A =PA 



dpi 

dt ' 



(25) 



where we have neglected the difference between p A and p c A in the first term on the right- 
hand-side. When this is nonzero it clearly makes no sense to define the time evolution in 
terms of the Schrodinger evolution, i-§i\ip) = H\ip), for a pure state: p A = \ip) (ip\. 



2.6) Recursiveness 

Notice that the definitions of the effective hamiltonian, V , and of the diffuse scattering 
term, dpi/dt, are recursive, in the following sense. Suppose that sector A in the previous 
discussion were itself to be divided into independent subsectors, A' and B' ', and only 
observables acting in subsector A' were measured. Suppose also that the initial state did 
not involve any correlations between these two subsectors: g A = q a > ® g B > . Then we may 
further reduce the density matrix to act only within this subsector: 



pA'{t) 



Tr 

B' 



p A {t) 



= Tr 

s'Ub 



Pit) 



(26) 



Then we define the coherent and diffuse part of the evolution in sector A' in such a 
way as to ensure that the coherent part takes the same form regardless of whether the 
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trace over sectors B and B' are performed separately, or all at once. That is, with the 
definitions 



U, 



Tr 

B 



QbU 



and U j 



Tr 

B> 



Qb' U b\ = Tr Uq b <S>Q B ')u], (27) 

J bus' L J 



and 



we have Tr, 



AU B = U-U B and AU B . = U B -U B ., 



(28) 



g B AU B = Tr B / £> S /A£/" B /J = 0, and so the expectation of any observable in 



only sector A' may be written: 

{0 Af )(t) = Tr[p(t)0 A , 

U B (t,t')g A Ul(t,t')G A/ 
U B ,{t,t')Q Al U* B ,{t,t')0 A , 



Tr 

A 



+ Tr 



AU B (t,t')p(t')AU* B (t,t')0 



A' 



= Tr 

A' 



+ Tr 



AU B ,(t,t')Q B ,AU*,(t,t')0 A , 



(29) 



+ Tr 



AU B (t,t')p(t')AU*(t,t')0 A , 



A similar expression holds for 



dp A > 

dt 



v B >(t)p A ,(t)-pAt)v B ,(t) 



dp d B , , dp 



dt dt 



(30) 



where: 



p%, =Tr AU B ,(t,t'W)AU* B ,(t,t') 

b' . 

p d B = Tr \AU B {t,t')p{t')AU%{t,t!)\ 

b'UbL ' J 



(31) 



The effective hamiltonians, V B (t) and V B >(t), are respectively defined, as usual, in 
terms of U B {t,t') and U B >(t,t'), using eq. (11). In particular, using the notation (• • ■) B , = 
l^s'tC' ' ')Qb'\i we have the following very useful perturbative expression: 



V B ,(t) 



((V(t)) B ) B , - i dr{ [({V(t)V(r)) B - (V(t)) B (V(r)) B ) B ,] 
+ [((V(t)) B (V(t)) b ) b , - ((V(t)) B ) B , ((V(r)) B ) B ] }, 
((V(t)) B ) B , - i fdr{{{V{t)V{r)) B ) Bl - ((V(t)) B ) B , ((V(r)) B ) B , }. 



(32) 
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The recursive nature of these definitions is a very attractive feature. This is because it 
lends itself to a renormalization-group-like analysis of the effects of a medium on particle 
propagation, in which the effects of fluctuations on successively larger distance scales are 
separately integrated out. 

2. 7) Neutrino Interactions 

Our later applications use the formalism just presented for the special case of neutrinos 
interacting with matter through the weak interactions. (Magnetic moment interactions are 
briefly considered in Section 7.) We therefore pause here to gather the relevant expressions 
for the interaction hamiltonian. None of the details of the nature of the medium are needed 
at this point, although we suppose for simplicity that it does not include a significant 
component of neutrinos themselves. 3 This permits a clean separation between the neutrino 
states and the states which are available to the medium. 

The neutrinos may be described by N v majorana neutrino fields, i>i,i = 1,...,N U , 
without loss of generality. In the absence of light sterile neutrinos this consists of the 
usual (N u = 3) neutrino eigenstates. The coupling between neutrinos and the medium is 
mediated by the weak interactions. To keep as broad as possible the applications of this 
section we choose: 

C = i V ilfJi ( 7i <$ + lR h a ij ) uj J£, (33) 

where 7i and 7h project onto left- and right-handed spinors; J£ are a set of hermitian 
operators involving the degrees of freedom of the medium, and gfj and hfj are corresponding 
N v x N v matrices of couplings. The reality of C implies gfj and hfj must all be hermitian. 4 

Since the most important applications are to the Standard Model (SM), possibly 
supplemented by various sterile neutrinos and/or neutrino masses, we record explicit ex- 
pressions for the quantities, and hfj in this case. The couplings to charged leptons, 
£ m , are given by: 

(J^) mn = zZ m7 ^(l±7 5 )C, 
with (g^ n ) ij = V2G F 

{g™)i3 = V2G F N tJ 4, 
and (/i± n ) lJ = 0, 

Although more than adequate for applications to the sun, this assumption can break down for su- 
pernovae or the early universe. 

4 Our conventions are tl =5(1+75), 7b=5(1— Ts) and V=iv^^°=—w^^o. 
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1 2 

2 + 4 



f34) 
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where V m i are the leptonic CKM matrix elements which arise in the charged-current cou- 
plings once neutrinos acquire masses, and Mij are the analogous matrices which can arise 
in the neutrino neutral-current couplings. s w denotes, as usual, the weak mixing angle. 
The corresponding couplings to hadrons arise through their quark content: 



(j£) a = iq a r(l±To)Qa, 
with ) y = y/2GMj (T 3a - Q a s 2 w ) , 

(gDii = V2G F M l3 {-Q a s 2 w ) , 
and (h±)ij = 0, 



(35) 



where T 3a and Q a are the third component of weak isospin and electric charge of the 
corresponding quark. 

For many practical applications the energies involved are sufficiently low that the 
relevant hadronic degrees of freedom are just protons and neutrons. In this case we may 
approximate eqs. (35) by the following effective macroscopic currents: 



(9 a + h (J©« + isDij (Jt)a « ^/=Mj [(1 - 4sl) iprfp - iu^n] 

+ (axial-current and higher-derivative terms). 

(36) 

We now use these expressions to compute the quantities V and in two different 
regimes. We first consider the case for which the fluctuations of interest occur on scales 
which are microscopic in comparison to those relevant to neutrino propagation. This is 
followed, in Section 4, by a consideration of macroscopic fluctuations, for which correlation 
lengths are much larger. 



3. Microscopic Fluctuations 



Consider first the case of fluctuations having microscopic characteristic correlation 
lengths. For neutrinos in the sun this includes the thermal fluctuations among the particles 
making up the solar interior, and so the expressions obtained in this section may in this 
case be tested against standard results. 

Our goal is to compute the quantities V and -^f- of Section 2. In order to apply the 
formalism we must first split the system into sectors, A and B. We choose A = v to consist 
of all of the states in the neutrino sector, while sector B = E ('environment') represents 
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the sector describing the other particles — electrons, nucleons etc. - 
neutrinos propagate. 5 

At first order in the interaction, eq. (33), we have AU = and: 



through which the 



Vi = (V) E = -Jd 3 x iv a „ (7,4. + lR h a ij ) Vj (JZ) E 



(37) 



To this order neutrino evolution is simply described by replacing the interaction current, 
J£, with its mean, (Jg) E = Tr B g E J% , evaluated using g E (= q b ), which describes the 
initial state of the medium through which the neutrinos pass. 

More can be said about the mean currents, j%(x) = (Jj£(x)) E , given more information 
about the state g E . With the sun in mind we take this to describe a mixture of nonrelativis- 
tic electrons, protons and neutrons which are mutually interacting dominantly through the 
electromagnetic and strong interactions. This permits the use of parity invariance to limit 
the form taken by the mean currents. Moreover, we also work in an electroweak basis, for 
which jVij = Sij and V m i = 5 m i for the three usual neutrinos, and j\fij = V m i = for any 
light sterile neutrinos. (The influence of neutrino masses is described in more detail once 
the trace over neutrino momentum states is performed in Section 5.) With these choices 
we obtain the usual estimate for the mean matter currents. We find (h%jJjf) = 0, and: 

(4 J a) E « ^§ {2** S je jt!(x) - S^(x) + 6^ (1 - Asl) [#(*) - #(*)]}, (38) 

where jf? = ze7 M e, j£ = ipj^p and j£ = in^^n are respectively the local electron, proton 
and neutron currents. This expression is obtained by averaging eqs. (34) and (35), using 
eq. (36). The axial- vector parts of the weak currents drop out of eq. (38) by virtue of parity 
invariance of the solar medium. Notice that, although the neutron current distribution is 
independent of the others, for slowly moving particles local electric neutrality implies 
j£(x) — j£(x), and so the terms in the last equation which are proportional to (1 — 4s^) 
cancel. For nonrelativistic particles there is also a further simplification, since we may 
neglect of all of the spatial components of the mean currents: j% xs n a 5£. 

At second order in the weak interactions two new things happen. First, V acquires 
a second-order correction, typically introducing to it an antihermitian part. Second, -4^- 



Our use of subscripts V and l E', in place of l A' and is meant to avoid confusion with the different 
choice for A and B which is made in the subsequent sections. 
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becomes nonzero, introducing decoherence into any propagating neutrino state. We find: 



V 2 = i [ dr d 3 x d 3 x' v ilpj (7^?. + lR h\ )) uj 
Jt' 



(39) 



and 



dpt 

dt 



drd 3 xd 3 x' jz/^ (7^ + 7^) i/j 
Pu Klx {l,9ki + lnh b M ) v[ (5J b x '5JZ) E + (*<-> r)}. 



(40) 



In these expressions p^(= p A ) denotes the neutrino density matrix at time t. A prime on 
any field denotes that it is evaluated at the spacetime point (x', r) — e.g. v' = z/(x', r) - 
while all unprimed fields are evaluated at (x, t). 

At this point we use the information that the scale of the fluctuations, £, are micro- 
scopic in comparison with the distances of interest for neutrino propagation. This means 
the correlations, (^5J^5J b x '') , can be written in the approximate form 

(6J£6J?) B MCg6 3 (x-x!), (41) 

where the coefficient functions, C£ b , are explicitly calculable given the state g E — a point 
to which we return below. 

Eq. (41) permits eqs. (39) and (40), for V2 and ' ^° ^ e wr itten as follows: 

V 2 = i [ dr d 3 x V ilt > {l.gfj + l R h%) Uj V' klx { lL g b kl + lR h h kl ) v' u (42) 
Jt' 



and 
V dt , 



2 



-2 f dr d 3 x C% 77 i7m (7,4 + lR h\ )) Vj Pv V' klx { lL g h kl + lR h b kl ) v[. (43) 
Jt' 



Notice that these interactions describe processes, such as vv — > vv or v — > vvv, in 
which neutrinos scatter from medium-dependent fluctuations. Similar interactions are fa- 
miliar for electromagnetic propagation through matter, where the analogs of eqs. (42) and 
(43) are quadratic in the electromagnetic field and so describe the scattering of electro- 
magnetic waves by microscopic fluctuations [14]. (Similar quadratic terms also arise for 
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neutrinos at second order in their charged-current interactions with the particles in the 
bath.) Since the matter-dependent effects depend differently on neutrino energies than do 
the same processes in vacuo they can, in principle, be separated from one another. 

The potential size of the matter-dependent correlation coefficients, C^ b A , may be esti- 
mated by computing them for thermal fluctuations in a system of nonrelativistic particles 
in local thermal equilibrium 6 and for which, for simplicity, we imagine there is a single 
conserved particle number, Af = f d 3 x J°, to whose current the neutrinos couple. 

For thermal fluctuations we compute all averages over E using the density matrix, g E , 
of the grand canonical ensemble: 

q e = Z- 1 e -(^-MAO/T (44) 

Z here is the standard normalization constant: Z = Tr E 6 -( h e-^)/t _ M ore generally, we 
consider media which are in local thermal equilibrium, and so for which g E has a similar 
form, but with mean thermodynamic properties which vary (over macroscopic distances) 
from place to place. The grand canonical ensemble is the one which is locally appropriate 
for this case [15], since the number of particles in any local region of the medium is not 
fixed. 

With these choices we may compute the local fluctuations of Af. Using the assump- 
tion that the constituents of the medium are nonrelativistic, we neglect all but the time 
component of the current: (J fI (x)) E = n(x) 5q. We find [16]: 



n 2 K T T + f £ (x) , (45) 



where T(x) is the local temperature, and k t (x) is the system's specific isothermal com- 
pressibility: k t = ^ {jj^j ' where p is the pressure. 

In eq. (45), fe(x) is a function whose scale of variation is the microscopic fluctuation 
length, £, and which satisfies the defining condition: J d 3 x fe(x) = 0. Because of these 
conditions, fi(x) can be neglected for macroscopic applications, such as when eq. (45) is 
used in eqs. (42) and (43) to describe neutrino evolution over scales which are much larger 
than i. 

Substituting eq. (45) into expressions (42) and (43), and using the results in eq. (10), 
reproduces the usual expressions [6] [7] [8] for neutrino scattering from a thermal ensemble 



6 This application presumes not being near a critical point for which thermal fluctuations need not be 
microscopic in size. 
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of particles. This can be seen by using the ideal-gas equation of state, p = nT, in which 
case k t = -£r£. It follows that the combination n 2 K T T — which through eq. (45) governs 
neutrino scattering — simply reduces to the particle density, n. This leads to a neutrino 
scattering rate, r ~ an ~ G 2 F mEn, which for solar neutrinos in the sun's centre (E ~ 1 
MeV and n c ~ 10 26 /cm 3 ) scattering from nucleons (m ~ 1 GeV) is negligibly small: 
r" 1 ~ 10 10 km. 

4. Macroscopic Density Fluctuations 

Our second application of the formalism of section 2 is to macroscopic variations in 
the mean currents, (J£) E , that arise in eq. (37). We do so partly because this source 
of fluctuations has until recently been ignored in the literature. More importantly, this 
type of fluctuation can produce effects which are much larger than those which arise 
microscopically. There are many situations in electromagnetism for which macroscopic 
fluctuations can furnish the dominant medium-dependent effects. A familiar example is 
furnished by the case of light propagating through a cloud. In this case the cloud is opaque 
because of density fluctuations on the scale of the water droplets which make up the cloud, 
rather than fluctuations on more microscopic scales. 7 Our purpose here, and in Section 4, 
is to analyze the analogue of such fluctuations for neutrinos. 

Eq. (45) implies that thermal fluctuations have negligible effects for neutrinos passing 
through the sun, so our starting point is the mean hamiltonian, eq. (37), which describes 
neutrino propagation after averaging over microscopic matter fluctuations: 

V E = - J d 3 x f a {x) iV lltl { lL g% + lR h a ij ) Vj , (46) 

where we write j£(x) = (J£(x)) E for the mean current. 

Now comes the main point. For any fixed current profile, j£(x), the propagation of 
any particular neutrino through the sun is perfectly well described by pure-state evolution 
using the mean hamiltonian given in eq. (37) (or, equivalently, eq. (46)). (This point is 
demonstrated in detail in Section 5.) It is, however, not in general true that successive 
neutrinos see the same profile, j%(x). On the contrary, successive neutrinos arriving at 
a detector may have been produced at different places within the sun and so can pass 

Of course, this analogy can be misleading if applied too literally to neutrino physics, since the absence 
of multiple scattering precludes neutrinos from 'refracting' from a large scale density fluctuation in the 
sun in the same way that light refracts through a water droplet. 
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through entirely different density profiles while en route to the earth. Alternatively, the 
density profile itself can change in the interval between the passage through the same region 
of different neutrinos. As a result, the neutrino flux to which a detector is exposed can 
be thought to have been processed through a constantly changing kaleidoscope of density 
profiles. 

We wish to adapt the formalism of Section 2 to describe the influence on neutrinos of 
this eternally varying current profile. We do so by modelling these density variations as 
being random in character. We therefore consider passing neutrinos through an ensemble of 
density profiles — whose properties are elaborated below — over which we must average 
to obtain the neutrino signal as seen by a detector on earth. Taking advantage of the 
recursive nature of the formalism of Section 2, we may simply take these formulae over in 
whole cloth, but with the mean hamiltonian of eq. (37) now interpreted as the microscopic 
hamiltonian, and with the averages over sector B = £ ('Ensemble') now interpreted as 
ensemble averages. 

The results are immediate. To first order in V E , the mean hamiltonian after the 
ensemble average now becomes: 

Vi = - I d 3 x (j£(x)) e iV^ ( 7l 0« + lR h%) v 3 . (47) 



Similarly, at second order we find: 



V 2 =i [ drd 3 xd 3 x' v ilfx { lL g^ + lR h a ij ) v 3 
Jv 



(48) 



£ ' 



—j =- J t/ dr d * x d3x ' fan {1l9?j + IrKj) Vj 



and 

'dpi? 

2 Jv k (49) 

Pu Pfc7A (7*0*1 + lnh b kl ) v[ (8fi'8ft) e + (t~ r)}, 

where 6j£(x) = j£(x) - {j£(x)) e . 

What remains is to estimate the ensemble averages, {j£(x)) £ and (8j£ (x)Sj£(x') /f . , 
which appear in these expressions. A key difference between these averages and those 
considered previously is that we may no longer assume the currents to be delta-correlated, 
as in eq. (41). 
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4.1) The Ensemble Properties 



The precise nature of these ensemble averages depends on the kinds of physics which 
is responsible for the varying currents that successive neutrinos see. It is useful to have a 
systematic framework within which to couch our later models of these fluctuations. This 
section outlines such a framework. 

Suppose, then, that the currents j£(x) are expanded in terms of a complete set of 
orthonormal functions, (f) N (x), as follows: 



i + ^2c aN 4> N (x) 



N 



(50) 



We take the coefficients, C aN , to be random variables having vanishing mean, which are 
uncorrelated for different modes: 



(Ccln) £ — and (C aN Ci, M ) ^ — C a &jv &nm- 

(51) 

This implies the following for the density distributions themselves: 

(fa{x)) £ 

{^{x)8f b {x')) e 

Notice that the completeness of the basis functions implies that the currents, j^(x), become 
microscopically correlated, (5j£(x) Sj^(x')) £ oc 5(x — x'), if C 2 ahN should be the same for 
all N. 

Clearly it cannot be true that the variables, C aN , are uncorrelated for all choices of 
basis functions. Different physical origins for the underlying randomness can lead to a 
different choice for the preferred, uncorrelated, basis. In what follows we use the following 
two models for the density fluctuations in the sun. 

• Locally- Varying Density Fluctuations: 

As our first model of solar fluctuations we picture j% (x) to be varying randomly from 
place to place. Motivated by the picture of the solar medium consisting of turbulent regions 
of fluid we imagine dividing the sun into cells, labelled by the index N, whose volume we 



fa ( x ) 3b( x ' ) Yl C "bN 0n (x) (f) N (x). 



(52) 



N 
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denote by V N . We permit this volume to vary slowly as one moves around the sun. We 
then choose the basis functions to be: 



4> N {x) 



V N 2 if x lies within cell N, 
otherwise. 



(53) 



With this choice we find the correlations: 




if x and x' both lie within cell N, 
otherwise. 



(54) 



We introduce here the dimensionless quantity e a b N by: C 2 bN = e^ bN V N , to remove the. 
dependence on the cell volume, V N , which enters due to the normalization condition for 
the basis functions, <f> N (x). 

How big might e a bN and V N reasonably be expected to be? For solar applications the 
convective zone is known to contain density variations on many scales [17]. Granules on the 
solar surface are ~ 100 km across. Giant convection cells are believed to have dimensions 
which are comparable to the depth of the convective zone itself: ~2x 10 5 km. Of more 
interest for neutrino propagation are the scales at the depths where neutrinos are produced, 
and where they resonate. Unfortunately, both of these regions lie within the radiative 
zone, where intuition based on the convective zone is unlikely to apply. Our analysis 
in subsequent sections of how these fluctuations modify MSW oscillations indicates that 
e 2 £ > 100 m is the range which is likely to have phenomenologically interesting implications 
for neutrino oscillations, where £ is the length of a typical cell in the direction of neutrino 
motion. 

This type of random model is very much in the spirit of refs. [10] and [12], for which 
the electron density is modelled as a random variable that is delta-correlated in space. In 
fact, eq. (54) directly reduces to the ensemble used in ref. [10] in the limit of negligible 
correlation length (taken in ref. [10] to be £ = 10 km), and when e 2 £ is taken to be constant. 
Ref. [12] makes a slightly different choice, ensuring a small correlation length by continually 
adjusting £ to be a tenth of the neutrino matter oscillation length, as this varies throughout 
the sun. Differences between our results and those of ref. [12] do arise for some regimes, 
which we believe to be due to this difference in treatment of the correlation length. 

Besides not assuming negligibly small correlation lengths, a more important difference 
between eq. (54) and refs. [10] and [12] is that we may take the fluctuations to vary 
differently as a function of position and time, as may be appropriate for some kinds of 
solar physics. Our next example presents an illustration of such a case. 
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• Helio 'seismic Waves: Oscillatory Normal Modes: 

Our second simple model of fluctuations in the sun is meant to model helioseismic p- 
waves. As a source of fluctuations through which neutrinos propagate, helioseismic waves 
have the great advantage of actually being known to exist. Furthermore, a fair amount is 
known about the spectrum and amplitude of these waves [18]. 

In order to strip away as much extraneous detail as possible, we start here with a 
simple wave within a rectangular geometry. (We report elsewhere on the results of more 
detailed modelling of neutrino interactions with helioseismic waves [19].) That is, we choose 
our neutrinos to be moving up the positive z axis, through a medium whose length in the 
z direction is L = 2R Q . A basis of modes which vanishes at the boundaries of this volume 
is given by: 



where we take the period and wavelength to be related in terms of the speed of sound, c s . 
To start with we take c s to be a constant, but we also present some results with c s = c s (z) 
chosen to more accurately mimic the properties of the sun. For constant c 8 , momentum is 
conserved in the z direction, and £± are determined by the boundary conditions to be a 



Finally, N± = N±(r±), denotes the dependence of the modes on the two other di- 
rections, r± = (x,y), transverse to the direction of neutrino motion. For a ^-dependent 
speed of sound, c s (z), we take N± to also be plane waves, as in eq. (55), labelled by the 
conserved transverse momentum, = k x e x + k y e y . For slowly varying wave amplitudes 
we then take the wave number in the ^-direction to be given by 



with u = 2tv/t±. In our subsequent numerical applications we take the wave amplitude to 
vanish for those z for which the resulting k z is imaginary. 

There are several types of physics which might be expected to produce a normal modes 
with a randomly varying amplitude. First, even if the sun were to be oscillating with a 
single mode, the amplitude of this mode as seen by successive neutrinos would differ. This 
is because, although any one neutrino sees an essentially static density profile, (since the 
time scale for neutrinos to entirely escape the sun is quite short — £ < L ~ several seconds 




(55) 



positive integer: £ + = L/(n — |) and £- = L/n. 




(56) 
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- compared to typical wave periods — r ~ several minutes) successive neutrinos can 
catch the wave at differing points in its cycle. Neutrinos passing through at random times 
would therefore see a wave with a randomly varying amplitude. 

Of course, the real sun does not simply ring with constant amplitude because various 
(poorly understood) processes permit energy to be transferred into and out of the various 
normal modes. This leads to additional randomness to the mode amplitudes, as seen from 
the neutrino's perspective. 

5. Tracing Over the Neutrino Momenta 

We now fill in the neutrino part of the picture, and compute how the fluctuations 
considered in the preceding Sections can enter into single-particle neutrino evolution. To 
this end we again apply the results of Section 2 to trace over the momentum/position part 
of the single-particle neutrino Hilbert space, with the goal of deriving the explicit form 
for eq. (10) acting in neutrino flavour space, with which we can analyze resonant neutrino 
oscillations. 

Two issues must be borne in mind when applying the results of Section 2 in this 
way. First, since we wish to keep track of the second-order effects, described above, due 
to matter fluctuations, we must use the recursive form for the effective description which 
was given in section 2.5. We therefore divide the neutrino sector itself into two subsectors, 
A' and B 1 ', with B' consisting of the span of all of the momentum states of the neutrino 
sector, whilst A' comprises the sector labelled by neutrino spins and flavours. To avoid 
confusion with our earlier choices for A and B, we introduce the new notation A' = F 
('flavour') and B' = P ('position') for this part of the analysis. 

Second, since all practical measurements of neutrino flavour also involve a position 
measurement — i.e. neutrino x is measured to be of flavour y when it arrived at point z 

- we must also remember to adopt the formulation of Section 2.2, in which we perform a 
partial measurement in sector B' = P. This innocuous point has important implications 
for the form of the fluctuation terms in the evolution equations for the reduced density 
matrix in flavour /spin space. 

To proceed we must choose the density matrix which describes the initial neutrino 
state. Assuming that the neutrino flavour /spin sector is initially uncorrelated with the 
neutrino momentum, g v = g F <8> g P , we must choose an explicit form for g P , in order to 
evolve the spin/flavour state, p F , forward in time. For a single neutrino, we would take 
this to be a pure, single-particle state, g P = describing an outgoing spherical wave 

packet which starts at t = t' = at the nucleus whose fusion produced the neutrino. The 
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spatial width, £, of this packet we imagine to be of negligible, microscopic, dimensions. 
Since applications to neutrino oscillations involve observations of this wave a very long way 
away from its centre, it suffices for our purposes to dispense with the spherical geometry 
and treat ipk,z (p) as a plane wave packet, starting at z = zq at t = 0, and travelling along 
the z axis with average momentum, k ~ MeV: 

(P\lpk,z ) = tpk{p) = 

iV'fc.zo) so defined is continuum normalized in the x and y directions. For solar neutrinos 
this pure state must be averaged over the initial distribution for producing such a neutrino 
within the sun. 

Next we must define the partial neutrino position measurement, which is meant to 
express the fact that we know where neutrino measurements are performed: the Earth. 
We therefore choose observables of the form (13) (O = O f <g> P ) with 

P (t m ) = \r;t m )(r;t m \, (58) 

corresponding to neutrino detection at the point, r, at a measurement time t = t m . For a 
long-term exposure to a constant flux (such as for solar neutrinos) we integrate over the 
appropriate range for t m . 

5.1) First-Order Effects 

With these choices we may now evaluate the quantity V\. (Second order effects due to 
the averaging over neutrino momenta are explored in the next section.) This will reproduce 
the usual MSW hamiltonian. The proper description of neutrino scattering, including the 
effects of matter fluctuations, is therefore found by tracing eqs. (47), (48) and (49) (or, for 
microscopic fluctuations, eqs. (37), (39) and (40)) over the neutrino momentum sector. 

This trace is straightforward to perform, subject to two important approximations. 

• Negligible Neutrino Masses: The first of these is the assumption, previously encountered 
in Section 3, that all neutrino masses may be neglected when evaluating V\. This is a good 
approximation for the masses and mixings which are relevant for solar neutrino oscillations. 
(Of course, neutrino masses do play an important role once V\ is used to evolve neutrino 
states forward in time.) 



2£ 



< 2 (Pz-k) 2 -ip z z 



${PX) S(j>y). 



(57) 
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• Slowly- Varying Density Profile: Secondly, j£(x) is assumed to vary negligibly over dis- 
tances comparable to the packet width, £, and to the neutrino wavelength, = 2n/k. 
This approximation implies the only significant scattering from macroscopic fluctuations 
is in the forward direction. For applications to solar neutrinos j£ varies macroscopically 
while £ and A„ are of atomic dimensions, so this last approximation also holds extremely 
well. 

We find the following effective hamiltonian, to first order in V: 

\Vi{r,t,t m )}iX;ja = -i J d 3 xd 3 q ^ )Zo (r, t m )^ Mo (q) j£(x,t) 

(r,X,i;t m \ u-/^ L g a + l R h a )v \q,a,j) (59) 
«n a (r,*) [M% 9(-a)-M?* 6(a)] 8 Xa . 

This result acts trivially on the spin labels, A and a, involving only the step function, 
9(a) = j[l + sig n(J ]? which projects onto left-handed (LH: a = — |) and right handed 
(RH: a = +|) states. We label the spin space using the projection of the spin in the 
z (or propagation) direction. For massive neutrinos this choice is made in the neutrino 
rest frame, while for massless neutrinos it applies in any Lorentz frame. Since A = — | 
corresponds to a left-handed state, we see that the within the Standard Model the states 
for which A = -h| are antineutrinos. For N u neutrino species the N v -hy-N v matrices 
represent the action of V\ on the flavour indices, % and j. They are given explicitly in 
terms of the coupling matrices by Ai^ = gfj — = gf- — hi* . 

Finally, the quantity n a (r, t) in this equation denotes the following: 

—0 —2 

n a (r, t) = j a [r± ,z + vt,t}- j a [r± ,z + vt,t], (60) 

where r± = (r x ,r y ) is the measurement position transverse to the neutrino propagation 
direction, and v denotes the speed, v = k/E^ ~ 1, associated with the central momentum, 
k, of the wave packet. Recall that zq denotes the point of origin of the neutrino, which is 
to be averaged at the end of the calculation. 

Using eq. (38) for the mean currents, j 1 ^, together with the nonrelativistic approxi- 
mation, which permits the neglect of f a — certainly good for the electrons and nucleons 
within the sun — we see that eq. (59), is recognizable as the standard MSW starting 
point for analyzing resonant neutrino mixing in matter. With this encouragement, we now 
proceed to compute the second-order contributions to neutrino evolution. 
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5.2) Second-Order Contributions 



The second-order contributions come in two kinds, as can be seen from eqs. (30) and 
(32) of Section 2.5. First and foremost, there are the matter fluctuations — i.e. eqs. (48) 
and (49) due to macroscopic fluctuations in the ensemble. But, a priori, there could also 
be diffuse contributions due to the neutrino momentum average itself, as was first discussed 
in ref. [5] . 

In this section we apply the general treatment of Section 2 to both types of fluctuations. 
For matter fluctuations we do so using the two ensemble models which were introduced 
in Section 4.1. The effective hamiltonian we obtain in this way will turn out to have 
interesting implications for MSW oscillations in the next section. By contrast, we find no 
phenomenologically interesting effects for solar neutrinos due to fluctuations which arise 
due to integrating out the neutrino momenta. Since this conclusions differs somewhat from 
that of ref. [5], we reproduce his results in our formalism, and show why our conclusions 
differ. 

We start with the formalism of Section 2, with the following three sectors: (i) A' = F 
for neutrino spins and flavour; (ii) B' = P for neutrino momenta and position; and (Hi) 
B = S (or E) for the matter degrees of freedom. Using the same approximations as were 
used to obtain eq. (59) for V\, we find the second-order contribution to to have the 
following form, regardless of the source of fluctuations: 



(V 2 ) iX , j(7 » -iA ab (r, t, t m ) (M a M b )ij 9(-a) + (M a *M b *) l3 9(a) 



'dpi 
, dt 



| » 2A ab (r, t, t m ) [(AT p F M b )ij 9(-a) + (M a * p F M h *) tJ 9(a) 

i\;ja 



(61) 



(Recall here r represents the position where the neutrino is detected at time t = t m .) 
The two kinds of fluctuations discussed above differ only in their predictions for the key 
coefficient, A a b(^, t, t m ). We now give expressions for this quantity for each of the two 
cases. 

• Fluctuations due to Tracing out Neutrino Momenta: 

Before computing A a b for matter fluctuations, we briefly pause to discuss the fluctu- 
ations which arise on integrating out the neutrino momentum sector. We do so partly to 
make explicit the contact with ref. [5] . We also do so partly because such fluctuations can 
arise and may be important in some circumstances. We argue here why solar neutrinos 
are unlikely to be one such case. 



28 



A direct application of the formulae of Section 2 to this type of second-order contri- 
butions shows them to vanish, within the approximations outlined above. This is because 
we have chosen to measure the neutrino position arbitrarily accurately, at precisely one 
point, (r, t m ). It is therefore instructive to integrate the observable, O p of eq. (58), over a 
finite detector volume, V. That is, we now replace eq. (58) by: 

O p = [ d 3 r \r,t m ){r,t m \. (62) 
Jx> 

With these choices, formulae (30) and (32), when applied to the neutrino-sector trace, 
give eqs. (61), with 

A a b(r,t,t m ) px dr (5n a [r ± ,z + vt,t]5n b [r±,z +VT,T]} p . (63) 
Jt> 

Here n a is as defined in eq. (60), and the average is over the detector volume transverse to 
the neutrino momentum. For any quantity, A(r,t), this average is defined by: 

(A(v,t)) p =^- [ d 2 r ± A(r ± ,r z ,t), (64) 
D± Jt> ± 

where V± denotes the area which the detector presents transverse to the neutrino beam. 
As usual, 5n a , denotes the deviation of n a from this transverse mean: 5n a [r, t] = n a [r, t] — 
{n a [r,t]) p . The key point here is that this deviation vanishes, Sn a = 0, in the limit that 
the transverse detector size, V±, is much smaller than the scales over which n a varies 
appreciably. 

This result makes sense physically. In the absence of matter fluctuations, neutrino 
evolution in the presence of a fixed density profile can be computed as an exercise in 
scattering from a fixed potential. Scattering only arises from variations, Sn, in the density 
from its spatial average. The main point is that the interference term between the scattered 
and initial waves in this problem is proportional to J v d 2 r± Sn(r±), and so vanishes only 
if the transverse area of the detector is sufficiently large on the scales over which 5n varies. 
It follows that the scattering is incoherent only for such large detectors. 

The relation of this result with that of ref. [5] is now clear. In this reference, the 
reduced density matrix for neutrino flavours is defined by completely tracing over all neu- 
trino momenta, without taking into account the position measurement, eq. (58). This is 
equivalent to taking the detector volume to fill all space, and our eqs. (63) and (64) indeed 
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reduce to ref. [5]'s in this limit. But working in the limit of an extremely large detector 
misses the important suppression of these effects by the detector size. 

Since we find in later chapters that significant neutrino effects require fluctuations on 
scales of hundreds of metres and up, we are led to conclude that this kind of neutrino 
incoherence likely plays no role for solar neutrinos. 

• Matter Fluctuations: 

For matter fluctuations, the second-order contribution to the neutrino evolution equa- 
tion is given by eq. (61), with 



A a b(r,t,t m ) = / dr (5n a [r ± ,ZQ + vt,t]5n b [r ± ,z + VT,T]) t 
Jt> 



(65) 



As before n a denotes the difference j® — j* (or simply the density, j®, in the nonrelativistic 
limit), while 5n a is the difference between n a , and its ensemble average, n a = (n a ) £ . 

Eq. (65) may be explicitly computed within the two models of fluctuations which were 
introduced in Section 4. For the case of locally- varying density fluctuations, eq. (53), we 
have: 

A a b(r,t,t m ) = e 2 ahN n a [r±,zo + vt,t] / n b [r±, z + vt, t] dr 

./cell N 

~ e 2 abN l N n a [r ± ,z + vt,t] n b [r ± ,z + vt,t] ^ 
(locally- varying density fluctuations). 

Here l N is the length of cell iV along the neutrino line of flight, and iV labels the specific 
cell which contains the point (r±, zq + vt,t). The approximate equality in the second line 
is derived under the assumption that the mean current, n a , does not vary appreciably over 
the size of this cell. For the approximate exponential density profile [2]: 

n a (z) = K) c e~ z /\ (67) 

which we use for electrons in the sun, the neglect of the variation of n a over a cell requires 
I < h = Rq/10.5 = 6.6 x 10 4 km. (Notice this is a much weaker condition than requiring 
I to be much smaller than neutrino propagation scales.) For electrons, the central density 
we use is (n a ) c = 1.5 x 10 26 /cm 3 . 

Similarly, A ab (r, t,t m ) may also be evaluated for the case of an oscillatory density 
profile. For solar applications, since r± is of order several minutes and (the light-travel 



30 



time across) £± < R Q is of order a few seconds, it is sufficient to neglect powers of £±/t± 
and Rq/t±. Again using the exponential profile, eq. (67), for n z (x), we find: 

•Aa b (r±,t,t m ) = (n a ) c (n b ) c e±(r ± ) F±[z + vt,t,t'], 

(68) 

(oscillatory density fluctuations) 

where e± = C±N± defines the dimensionless size of the fluctuation, and the function F is 
defined by 

rz +t 

F ± (zo + vt,t,t') = e-^ +vt ^ h f ± (z + vt) / dxe- x ' h f±{x) (69) 

with f-(z) = sin {^f^j and f+(z) = cos {^ff^j- The integral is elementary and is given 
(writing a = 1/h and b = 2tt/£) by: 

F_(* , t-t') = e ~ a{Z ° +t) s ™K z ° + *) L-a(zo+t') \ bcosb ( Zo + t ') + a smb(z +t>) 

a z + b z L L J 

F + (z + t, t, t') = e ~ a{Z0+ l 2 C ° S p Z ° + t] {e-^« +t ') [a cos b(z + t') + b sin 6(^ + 

(70) 

To summarize, fluctuations can indeed influence the neutrino evolution. Their effects 
are quantified by equations (61) and (65), which are the main results of this section. Their 
implications for resonant MSW oscillations can be sizable, as is now explored in more 
detail. 



6. Applications to MSW Oscillations 



In this section we evolve the neutrino density matrix to second order in G F . With 
solar neutrinos in mind we follow the usual practice and suppose the initial neutrino spin 
to be purely left-handed, and focus on the evolution in flavour space. The plan is to 
use eqs. (59) and (61) to evaluate the right-hand-side of eqs. (10) and (12), and then to 
integrate the result to determine the electron-neutrino survival probability, P e (t) = Pee(t)- 
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6.1) The Evolution Equations 



For simplicity we specialize also to the case of two neutrino flavours, whose electroweak 
eigenstates we denote e and fx, although we might equally well imagine mixing the e-type 
and r-type neutrinos. Eq. (23) then takes the following form: 



dp 
dt 



Vo + V u p 
Vo + V u p 



2G 2 F A ab (g a g b P + P9 a g b -2g a pg b ), 



a,b=e,n 



2 G 2 F A ee [(g e ) 2 p + p(ge) 2 -2g e pg< 



(71) 



where p = p F is the neutrino-sector density matrix in flavour space, and 



V = 



k 2 + rrvm 



1/2 



k + 



V x {t) = V2G 1 



2k 

g e n e (t) + g n n n (t) 



■ + 



(72) 



Here m, g e and g n are 2x2 matrices which represent the left-handed-neutrino mass matrix, 
and the neutrino charged- and neutral-current coupling matrices. In an electroweak basis 
these are given explicitly by: 



m = 



m e 
m 



m 
m 



e fi 



g" = 



1 o 
o o ' ' 



and 



9 = 







-■ 



(73) 



In general m may be a generic symmetric complex matrix, although in the absence of CP 
violation it may be chosen to be real. n e (t) and n n (t) are the spatially-averaged electron 
and neutron currents, as defined in eq. (60). 

The second-order contribution to eq. (71) is given by the three terms proportional to 
A a b(t), which is defined for a, b = e, n by eq. (65). The terms of the form g 2 p and pg 2 are 
the contributions due to V2, while the gpg term comes from Notice that because g n 
is proportional to the unit matrix, all but the term involving A ee — which we henceforth 
denote simply by A — give zero in the sum in the first of eqs. (71). As a result, it is only 
fluctuations in the electron density profile which are relevant for neutrino evolution in the 
sun. 

In order to integrate this equation it is useful to expand all matrices in terms of the 
unit and Pauli matrices, {I, r}. We have 



P = Po + P ' T an d 
g e =l(l + T 3 ), and 



V = M + M ■ r, 



(74) 
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with 



Pee = P0 + P3, 



Pw = PO ~ P3, 



Pe M = Pi - IP2, 



(75) 



and 



M = k + 



M, = — 



\m eil \ 2 +\ (|m ee | 2 + \m^\ 2 ) 
2k 

Im (m* ee m etl +m* etl m^) 



2k 



Mi 
M 3 



Re (m* e m eM + m*^m MM ) 



2A; 

I |2 I I 



(76) 



Of these components, Mo plays no role in the evolution of p since it drops out of the right- 
hand-side of eq. (71). Similarly, since the trace of the right-hand side of eq. (71) vanishes 
identically, it follows that the coefficient po is independent of time: = 0. 

With these definitions the evolution equation, (71), for the remaining three compo- 
nents of p may be written: 8 




-2a 
2(M 3 + b) 
-2M 2 



-2(M 3 + b) 2M 2 
-2a -2Mi 
2Mi 



(Pl\ 




P2 


= n\ P 2 




\P3 



(77) 



where 



a(t) = G 2 F A{t) and b(t) 



G F n e (t) 
V2 ' 



(78) 



This is the form for the evolution which are integrated in subsequent sections. 

For the purposes of exploring the implications of electron density fluctuations it suffices 
to restrict our attention to CP-conserving neutrino physics, for which the components of 
Vo simplify somewhat because the neutrino mass matrix, m^, may be chosen to be real. 
For this case Mi through M3 can be expressed in terms of the heavy and light neutrino 
mass eigenvalues, mh and mi, by: 



bra 2 sin 26 v , , n , , 5m 2 cos 29 v 
= 4k M2 = °- M3 = — 

where dm 2 = m\ — m 2 and 9 V is the vacuum mixing angle, for which 

v e cos6 v — sin6*v, 



(79) 



v e sin 9 v + cos9 v . 



(80) 



If A(t) is assumed to be a constant times n 2 (t), and in the absence of CP violation, then this equation 
agrees with that used in ref. [10]. It also agrees with ref. [12] if the correlation length is adjusted as explained 
in Section 4.1. 
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We now turn to integrating eq. (77). Before turning to a numerical solution using an 
exponentially falling electron density profile [2], we first pause for the instructive exercise 
of solving this equation in the adiabatic and Parke limits. 



6.2) Adiabatic Evolution and the Generalized Parke Formula 

If a and b are slowly varying, then it is straightforward to analytically integrate eq. (77) 
simply by performing a (t-dependent) rotation which diagonalizes the matrix 7i. Once 
these adiabatic solutions are obtained, then Parke's more general expression may then be 
derived by starting with these states as bases and computing the transition probability as 
the neutrinos pass through the resonance. 

The adiabatic result is: 




/><*) 



dx 



\ 









(:) 


/ 





(81) 



where Xi(t) are the three time-dependent eigenvalues of 7i, and R(t) is the matrix for 
which Kl(t)H(t) R(t) = diag (Xi(t), A 2 (£), A 3 (£)). 

Keeping in mind that a arises at second order in G F and so is smaller than all of 
the other elements, we may solve for the Aj and R perturbatively in a. The resulting 
eigenvalues are, to linear order in a: 



Aq = — 270, A± = ±2m — 27, 



(82) 



with 



K 



^M? + (M 3 + b) 2 , 70 



aM 2 



q[M 2 + 2(M 3 + 6) 2 ] 
7 = nr „2 rm—, r^n • ( 83 ) 



M 1 2 + (M 3 + 6) 2 ' ' 2[M 2 + (M 3 +6) 2 
To the lowest (zeroeth) order in a the corresponding matrix, R, is: 

i i_ n i i_ 

V2 V2 — u ^ ^ 

\-cos26 m -^sin20 m -^sin2(9 m y 



Mi 



M 3 +b 



/ 



where the last equality defines the usual matter mixing angle, 9 m . 



(84) 
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Using eqs. (82), (83) and (84) in eq. (81) then gives the adiabatic prediction for the 
electron-neutrino survival probability, given the initial condition p(t') = diag (1,0): 



p e {t) = p ee {t)=\ (l + e 2 J> ^ cos20 m (O cos2# m (t) 

" 1 (85) 



+e 2 I t ,-r(^ dx S i n 2# m (t) sin20 m (£') cos 



K,(x)dx 



The oscillatory term, on the second line of eq. (85), averages to zero once we sum over a 
long enough measurement interval. 

For nonadiabatic evolution it is straightforward to use these adiabatic results to derive 
a generalization of Parke's formula. After averaging over the production and detection 
times we find: 

P e (t) =1+ Q-Pj) e- 2 I> o(x)dx cos20 m (O cos2# m (t), (86) 

in which P d = exp — ^ i^W^) * s tne 'j um P' probability as one passes through 

the resonance point. 

There are many reasons to become emotionally involved with eq. (86): 

• 1: In the absence of fluctuations, a — > 0, eqs. (85) and (86) reduce to the standard results 
for matter oscillations. 

• 2: When a is small, but not zero, its dominant influence is to damp the neutrino oscilla- 
tions, by introducing an imaginary part to the masses of the mass eigenstates in matter. 
Such an imaginary contribution might have been expected given the antihermitian form 
found in eq. (61) for V 2 - 

The resulting damped oscillations are similar to what arises when neutrinos decay, 
but with an important difference. The difference is the appearance of the term in the 
evolution equation, (10). As a result, there is no net loss of probability in the neutrino 
sector: Trp F (t) = 1 for all t. For neutrinos the damping is a reflection of the conversion 
of the incoming neutrino from a pure to a mixed state, due to the decoherence introduced 
by the fluctuations in the matter through which it passes. 

• 3: The relative size of the fluctuation parameter, a = G 2 F A, in comparison with the usual 
MSW effective hamiltonian, b = G F n/ \/2, is given (for locally- varying fluctuations) by 
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which gives a rough indication how large fluctuations must be in order to contribute 
significantly to neutrino evolution. 

Important effects arise even for smaller e 2 £ : however, in the case of resonant oscilla- 
tions. This is because even a small damping term can ruin the quality of the resonance. 
This is borne out by our numerical integrations, which show that the first modifications 
arise precisely at the resonant point. It can also be seen analytically from eq. (86), as we 
now argue. 

• 4 : The strength of the decohering effect is completely parameterized by the integral over 
70 which appears in eq. (86), after the oscillatory contributions are averaged out. This 
integral is elementary (for an exponentially falling electron density profile) when A(x) is 
a constant times n 2 (x), as is the case for the locally- varying density-fluctuation scenario, 
discussed earlier. In this case we have A = e 2 £n 2 and so 



1Z = exp 



= exp 



exp 
exp 



7o(x) dx 



, * dvt 

exp { -2h / 70 (n) — 
Jo n 



e 2 £h (5m 2 ) 2 sin 2 29 



In 



sin 20, 



sin 29' m 



4k 2 

+ cot 29 v ( arctan cot 29 v — arctan cot 29' 



Tve 2 £h(dm 2 ) 2 sm29 l 



4k 2 



■X 



e 2 £ 



dm 2 jk 



7kmy \\ eV7 MeV ) \ 0.1 



l + 0(sin2# v ) 



/sin 29 x 



(88) 



Here t ex is the time the neutrino exits the sun, and so after which n e (t) and 70 (£) both 
vanish. 9' m = 9 m (t') and n' = n e (t') are the matter mixing angle and electron density at 
the production point, deep within the solar interior. 

Eq. (88) shows that significant damping is possible for solar neutrinos when e 2 £ > a 
few kilometers. 

• 5: For the smallest amplitude fluctuations, the first place where the damping becomes 
noticeable is for resonant, adiabatic oscillations (see Figures 1 and 2). That is, for Pj ps 
and cos 29' m ~ — 1, and taking 9 m (t) = 9 V , eq. (86) becomes: 



-T e ~ 



1 _ g -2 J t * -yo(x)dx 



cos 29 x 



(89) 
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For a = (no fluctuations) this gives the usual suppression: P e (t) « sin 2 6v- For small 
damping this suppression is weakened. 

Since the success of MSW oscillations in explaining the solar-neutrino data uses this 
suppression to virtually remove the 7 Be neutrino line, a good measurement of the strength 
of this flux promises to give information about the strength of electron density fluctuations 
deep within the solar interior. 

• 6: For sufficiently large times, J^, , yo(x)dx ^> 1, eq. (86) has the universal prediction: 
P e — > |. (This limit is also seen in our numerical integrations, as well as in those of 
ref. [10]. We believe its absence in ref. [12] is due to the use there of a correlation length 
which follows the neutrino matter oscillation length.) This suggests a new solution for 
the solar neutrino problem: an approximately energy-independent suppression of the solar 
neutrino flux by a factor of 2 due to solar fluctuations, even for small neutrino mixing 
angles in vacuum. This type of solution will be disfavored if improvements in the data 
should confirm the present indications for an energy-dependent neutrino suppression (i.e. 
for which p — p neutrinos are untouched while 7 Be neutrinos are essentially removed). 

• 7: Resonance is defined by the condition M3 + 6 = 0. Interestingly, although this 
resonance condition minimizes 7, it actually maximizes 70 This is clearest when eq. (83) 
is written: 70 (t) = a(t) sin 2 29 m (t). The sharper the resonance, the sharper the peak there 
in 7 . In the limit of a sudden resonance we therefore expect 1Z, of eq. (88), to be controlled 
by the value of A evaluated at the resonance point. In this limit we expect the fluctuation 
size at the resonance point to be what determines the size of the damping contribution. 

6.3) Numerical Results 

We have numerically integrated eq. (77) to determine the electron-neutrino survival 
probability, P e (t) = Pee(t), using A(t) as predicted by both the locally- varying density- 
fluctuation and the oscillating-mode model of density fluctuations. For the oscillating-mode 
case we have examined both the case of constant sound speed, and a speed which varies 
as a function of z. 9 We use the exponential density profile, eq. (67). 

The present section is devoted to presenting the results of these calculations. We have 
compared these numerical results with the analytical result, eq. (86), which we find works 
extremely well for all of the parameters of interest. 

9 We use the form c(z)=c 1 +c 2 (jr)+c 3 (j : y +czt(+) 3 +cs(-|-) 4 , in units of the speed of light, with coeffi- 
cients ci=0. 00170, c 2 =0. 000581, c 3 =-0.0106, c 4 =0.0177 and c 5 =-0. 00929, as found by fitting to the profile 
given in ref. [18]. 
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We present our results through Figures 1 and 2. These present plots of the survival 
probability, P e (t), as a function of the neutrino energy normalized to its squared mass: 
E/8m 2 . We use sin26> v = 0.1. In Figure 1 we use the locally- varying density fluctuation 
model, with successive curves representing different values for the fluctuation strength, e 2 £. 
As expected based on eq. (86), sizable effects first appear by deteriorating the quality of 
the suppression at resonance. The first deviations appear when e 2 £ = 0.01 km, until the 
resonance is completely destroyed for e 2 £ = 100 km. For small E/Sm 2 and e 2 £ < 10 km, 
the survival probability approaches its asymptotic MSW form P e cos 2 26 v 1. For e 2 £ 
larger than this value, P e instead tends to the limit \ for small E/5m 2 . 

Figure 2 presents a similar plot for an oscillatory fluctuation, using the position- 
dependent speed of sound. The wave period is chosen to be 30 minutes, and its transverse 
momentum, k^, is chosen to be as large as is possible: k± ~ 2n/L, for L ~ R & . The 
wavenumber, k z , in the direction of neutrino motion is fixed as a function of z using the 
dispersion relation, eq. (56), with a realistic ^-dependent speed of sound. As was described 
in section (4.1), the wave's amplitude is set to zero when k z so determined is imaginary, 
corresponding to a damped wave. Our plots are made using an 'odd' mode, 

The figure displays the resulting survival probability as a function of E/5m 2 for various 
amplitudes, e. Sizable deviations start for e > 1%. The size of this deviation roughly agrees 
with what would be estimated using the locally- varying density fluctuation results of Fig. 1, 
with the correlation length, £ taken as the wavelength at the resonant point. This leads 
us to expect a larger effect in a more accurate simulation using g waves than for p-waves, 
since these have larger amplitudes and wavelengths in the resonance region. More realistic 
simulations to investigate these issues are under way [19]. 

7. Magnetic-Moment Couplings 

This section briefly applies the formalism of Section 2 to neutrinos which interact with 
magnetic fields through a magnetic-moment interaction: 

£ = ^cTwA F» v + ex. 

1 r n ( 90 ) 



z^7^(Re nij + i>y 5 Im 



± [lis- 



Since the Vi are majorana, the matrix fiij is (complex) antisymmetric. Fluctuating mag- 
netic fields have been considered in ref. [10], in the limit of negligibly short correlation 
length. 
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With this choice, the first-order contribution to the effective hamiltonian governing 
neutrino evolution in spin and flavour space as it moves along the z axis is: 



(Vi)i\-ja = Tij sign(cr) S-X,a, (91) 
with Tij representing the following combination: 

T ij = Re Hij (E x + By) + Im ^ (E y — B x ) 



Re 
Re 



Hij(E x -iE y )^ -Im fjiij{B x - iB y )^ . (g 2 ) 



Here Ek and B\ denote the (microscopic or ensemble) mean electric and magnetic fields, 
averaged over the matter sector: 

Ek = (Fok)g and B k = -tklm {Flm)s > (93) 

and S k = E k + iB k . 

Similarly, the second-order contributions to the evolution equation become: 



{V2)iX-ja = Sxa [ dr ({5F6F) £ ) 

Jt' V 7 V 

sign(A) sign(a) / dr ((S^pp-SJ 7 )^ 

Jt' v J * a ;j <j 



where 8Tij = Tij — Tij. 

Clearly these expressions can be used to extend the analysis of fluctuations in magnetic 
fields to systems for which the fluctuations are uncorrelated in momentum space. Based 
on our experience with the MSW oscillations, for reasonably small fluctuations we expect 
appreciable consequences dominantly in the presence of resonant oscillations. 

8. Conclusions 

In this paper we have obtained the following results: 

• 1: We set up a general formalism for describing neutrino propagation through fluctuating 
media. This formalism has the virtue that it is derived from first principles, and so there 
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are no hidden assumptions which limit its applications. As a result it can be applied to 
any source of fluctuations which may be of interest. 

• 2: We have applied this formalism to many of the sources of fluctuations for neutrinos 
propagating through the sun. When applied to microscopic, thermal fluctuations we repro- 
duce standard estimates, which give a negligible impact on neutrino propagation. When 
applied to larger-scale macroscopic density fluctuations the effects can be larger, typically 
becoming important once e 2 £ > 100 m or so. Here e ~ dn/n is the relative amplitude 
of the density fluctuation, and I is a measure of its correlation length in the direction 
of neutrino motion. If this varies from place to place, it is its value at the point where 
resonant oscillations occur that is most important. 

• 3: The neutrino evolution equation which we obtain — eq. (77) — agrees, when restricted 
to the domain of common validity, with those of previous workers [6] [7] [8] [5] [10] [12]. We 
find the fluctuations found in ref. [5] to only significantly decohere the neutrinos if their 
size is small compared to the size of the detector. Since the strength of their influence on 
neutrino propagation is itself proportional to the fluctuation size, such small fluctuations 
are likely to be negligible for solar neutrinos. 

• 4 : Two models of macroscopic solar fluctuations were developed. One of these, having 
cells of constant, fluctuating density, reduces to ref. [10], and is similar to ref. [12], in the 
limit of small cell size. The other is completely different, and models the oscillatory density 
variations which occur in solar acoustic p-waves. (A more realistic version of this model is 
currently under study.) For these oscillatory waves, taking a period of 30 minutes, we find 
appreciable neutrino effects for density fluctuations which are at least a percent in size. 
We trace this comparatively small effect to the relatively small wavelength at the neutrino 
resonance point, which comes about because, for p waves, the wavelength decreases with 
depth due to the increase of the speed of sound. We do not expect the same suppression 
to apply to g waves. 

• 5: We integrate the resulting neutrino evolution equations for the case of two neutrino 
flavours (with no CP violation) to obtain a generalized Parke's formula — eq. (86) - 
for the electron-neutrino survival probability. This formula reproduces well our numerical 
integrations. The decoherence due to the fluctuations enters the neutrino evolution like a 
friction term, causing the oscillations to damp. 

Two features emerge from the result. First, for small fluctuations deviations from 
the usual MSW survival probability first arise for adiabatic resonant transitions, for which 
the MSW suppression deteriorates because the fluctuations partially ruin the resonance. 
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Because MSW oscillations use this suppression to agree with the solar neutrino data by 
eliminating the 7 Be line, measurements of this line in neutrino detectors promises to shed 
light on the nature of density fluctuations deep within the sun. 

Second, for small E/5m 2 and for sufficiently large fluctuations, the survival probability 
falls to 0.5, independent of energy. This introduces a new fluctuation-driven mechanism 
for solving the solar neutrino puzzle, although its energy dependence is not favoured by 
current measurements. 

• 6: Finally, our formalism is used to derive the effective neutrino hamiltonian which is 
relevant to magnetic moment couplings to magnetic fields. 
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10. Figure Captions 



(1) The electron neutrino survival probability as a function of E/5m 2 assuming the locally- 
varying density fluctuations model, as described in the text. Each curve represents a 
different value for the parameter combination e 2 £, where e = dn/n is the fractional am- 
plitude of the fluctuation, and I is the fluctuation's correlation length in the direction 
of neutrino motion. 

(2) The electron neutrino survival probability as a function of E/dm 2 assuming the os- 
cillatory density fluctuations model, which is a crude model of a helioseismic p-wave. 
Each curve represents a different value for the fractional amplitude, e = dn/n. The 
figure assumes a wave period of 30 minutes, and uses the speed of sound as a func- 
tion of depth given in ref. [18]. The transverse wavenumber, kj_, is chosen as small 
as possible, corresponding to a wave which penetrates as deeply as possible into the 
solar interior. 
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Figure 1 
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